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Abstract 

Soliton interactions in systems modelled by coupled nonlinear Schrodinger (CNLS) equations and 
encountered in phenomena such as wave propagation in optical fibers and photorefractive media 
possess unusual features : shape changing intensity redistributions, amplitude dependent phase 
shifts and relative separation distances. We demonstrate these properties in the case of integrable 
2-CNLS equations. As a simple example, we consider the stationary two-soliton solution which is 
equivalent to the so-called partially coherent soliton (PCS) solution discussed much in the recent 
literature. 
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I. INTRODUCTION 



The study on the formation of optical solitons and their intriguing interaction properties 
is becoming one of the frontier areas of research in nonlinear dynamics due to their potential 
technological applications [1,2]. Indeed optical solitons are becoming desirable candidates in 
long distance optical communication systems, in optical devices and in optical computers. 
In a mathematical sense these solitons appear basically as solutions of integrable coupled 
nonlinear Schrodinger (CNLS) type equations. For example, the intense electromagnetic 
wave propagation in a birefringent fiber is governed by the following set of 2-CNLS equations 
for the envelopes q 1 and q 2 , which is in general nonintegrable, 

iqiz + qitt + 2/i(|gi| 2 + B\q 2 \ 2 )qi = 0, 

%q2z + qm + Ml^f + B\ qi \ 2 )q 2 = 0, (1) 

where z and t represent the normalized distance along the fiber and the retarded time re- 
spectively, n represents the strength of nonlinearity and B = ^^iSe * s ^ ne cross phase 
modulation coupling parameter (9: ellipticity angle). However this system becomes inte- 
grable for B — 1. The resulting set of equations 

iqiz + qui + 2M|<?i| 2 + l<?2| 2 )<?i = o, 

iqiz + q-m + 2/i(|<?i| 2 + l<?2| 2 )<?2 = 0, (2) 

is the celebrated Manakov equation[3]. In a recent work Radhakrishnan, Lakshmanan and 
Hietarinta [4] have revealed the fact that the soliton solutions of the integrable 2-CNLS 
(Manakov) equations undergo a fascinating shape-changing collision, resulting in a redistri- 
bution of intensity between the two solitons in the two modes, which is not observed in the 
scalar nonlinear Schrodinger (NLS) equation which exhibits only pure elastic collision with- 
out any redistribution of intensities of solitons. Consequently, Jakubowski, Steiglitz and 
Squier[5] have pointed out the possibility of using this phenomenon in constructing logic 
gates and in a very recent work [6] Steiglitz constructed such gates including the universal 
NAND gate, thereby showing the theoretical possibility of constructing all optical comput- 
ers without interconnecting discrete components in a homogeneous bulk nonlinear optical 
medium. Also Yang[7] has studied the effect of additional perturbations on these solitons 
using perturbation theory. Further such integrable CNLS equations arise in the context of 
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spatial solitons as well which are receiving renewed attention for their formation at very low 
optical powers in photorefractive medium[8]. The present authors have extended the results 
of 2-CNLS system to 3- and N-CNLS equations [9]. 

All the above investigations mostly concentrate on the effect of changes in the amplitude 
(polarization) and the consequent effect on the energy redistribution between the modes of 
the solitons. So far not much attention has been paid to the role of phases during optical 
soliton interaction. In this report, we point out the significance of amplitude dependent phase 
shift / relative seperation distance involved in the interaction process (sec. 2) responsible 
for the shape change of solitons during collision along with the changes in the amplitudes 
(polarization) of the modes in the Manakov system. It may be noted that such amplitude 
dependent phase shifts do not occur in the case of scalar NLS equation. As a simple example, 
we consider (sec. 3) the role of phase shifts / relative separation distances for stationary 2- 
soliton case and point out that this solution is nothing but the so-called stationary partially 
coherent soliton (PCS) in the recent literature [10]. 



II. SOLITON INTERACTION IN 2-CNLS SYSTEM 

To start with let us consider briefly the nature of one- and two-soliton solutions [4,11]. 
Multisoliton solutions of CNLS equations and their interactions will be considered else- 
where [12]. 



A. One-soliton solution 



The one-soliton solution to eq.(2) can be given in terms of three arbitrary complex pa- 
rameters af\af^ and k\ as[4,ll] 

\ / ( ~\\ \ 

(3) 
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I _|_ e m+vl+R : 
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cosh (r] 



m 1- 2", 



HT> ( 4 ) 



where rji = ki(t + ikiz), i — 1, where h L = km + iku, km and kn represent the real and 
imaginary parts of fcj. Here v /u(Ai,A 2 ) = , = ' Ql represents the unit polarization 

V|a«|2+|4 2) |2 

vector, ki R Aj, j=l,2 gives the amplitude of the jth mode and 2ku is the soliton velocity. 
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B. Two-soliton solution 



Introducing six complex parameters a 2 \ ct[ \ a^, k\ and k 2 , the two-soliton solu- 
tion can be given as [4,11] 



where 



Here 
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and/Cij = \ ■ J - , i,j = 1,2. 



The above two-soliton solution represents the interaction of two coupled one solitons. The 
scenario behind this interaction is that there is an intensity redistribution among the two 
modes of the two solitons along with an amplitude dependent phase shift and relative sep- 
aration distance[4,9,ll]. In order to understand the nature of the collisions we can consider 
the following cases for ku>k 2I : 

(a) ki R > 0, k 2R > (b) k m > 0, k 2R < (c) k m < 0,k 2R > (d) k m < 0,k 2R < 0. 
Similarly, one can consider four cases for ku < k 2I . In all these cases an asymptotic analysis 
(z — > ±00) reveals the following structures. 



1) Limits — > — oo 
(a) Solitonl : 



(b) Soliton2: 
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2) Limits — > +00 
(a) Solitonl : 
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(b) Soliton2: 
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Here the various quantities corresponding to the four cases for ku>k 2 i are given 
Case(a)A;i_R > 0,k 2R > 0: 



A\- 
A\+ 
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and 2+ = ^. Case(b)A; li? > 0,k 2R < 0: 



-(R i + J R 3 )/2 
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ft' = and <^' + = X> 3 = 1, 2. 

Case(c) fcifl < 0,A;2R > 0: 
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= 5l an d = *a=*i. 
In both the cases (b) and (c) I = j + (— j 
Case(d) k 1R < 0,k 2 R < 0: 
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In equations(6-21) the superscripts denote the solitons and the subscripts denote the modes. 
It is clear from the above expressions that in all the cases there exists a redistribution of 
intensity among the solitons. However it should be noticed that though there is an intensity 
redistribution among the solitons in two modes, the total intensity of individual soliton is 
conserved during collision process, that is, |Ap| 2 + | ^4_2 | 2 = |^4{ + | 2 + |^2 + | 2 = ~> 3 = 1>2, 
which is of course a consequence of the integrability of the Manakov model. For example, 
the amplitude change in the two modes of soliton 1 after interaction can be expressed by 
the following transformation, 

A\ + = TC lx A\- + TC l2 A\~, 

A\ + = TC 21 A\- + TC 22 A\-. (22) 
6 



Here T = T(A\~ ,A\~ ,A\~ , A\~) = (g) [l/Ka^* + af^X^ 
r , , 1-V2 
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, in which a 2 = (k 2 + k^Kh - k^a^ a ( 2 ] * + af ] a^*)} 1 ' 2 . The forms of 



a^a^ih-^ + a^a^ih + k*), C 12 = -a^af^ + A;*), 



C 21 = -a i 2 2) a { 2 ) *(k 2 + k*), C 22 = a^a^ih + k 2 ) + c4 2) c4 2) *(£;i - k 2 ). Note that C^'s are 
independent of a^'s and so of A\~ and A 2 ~. Similar relations for the soliton 2 hold good 
for A 2+ and A 2 ^ also. Then the ratios of the A?'s, i,j = 1,2, can be connected through 
linear fractional transformations(LFTs). For example, for soliton 1, 

A\ + C llP r + Ci 



'12 



A\ + C 21 p{ + C : 



(23) 



22 



A 1 - 

where pi=^r=, ensuring that for every transformation there exists an inverse transformation. 
This idea has been profitably used in constructing logic gates [5, 6]. In fact the LFT (23) is 
identical to the LFT given by equation(9) in ref. [5] under the change of notation p^—tpi, 
Pi^Pr with Cy's identified as the expressions given therein. 

Further, it is observed that the absolute value of the phase shift of the two solitons in all 
the above four cases is same and is given by 
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where 
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), hj — 1,2. and the absolute value of the change in 
relative separation distance tf 2 (position of S 2 (at z — > ±oo) - position of Si (at z — > ±oo)) 
is given by 
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(25) 



It is interesting to note that in the collision process the phase shift is not only dependent on 
the k/s, j — 1, 2 but also on the complex parameters a-^'s, i, j = 1, 2, and so on A^s. These 
two properties, that is, dependence of the change in the intensity profiles of the solitons in 
the two modes and of the phase shift of them during collision on the parameters a^'s make 



the collision properties novel, not seen in general other standard (1 + 1) dimensional soliton 
systems. 

Now looking at the dependence of the phase shift on ce^'s, we can consider two special 
cases. 

Case(a): : 



(2) (2) 



In this case (corresponding to parallel modes) the collision correspond to pure elastic collision 
(\Al + \ = \A 3 ~\, i,j = 1,2) and the phase shift is given by 



l$l = 



log 



\h-k 2 \ 2 



(26) 



Case(b): : 0$ =oo, a(' : a 2 ' = 0. 



(2) . J 2 ) 



This case corresponds to two orthogonal modes. Here the phase shift is given by 

\ki - k 2 \ 



log 



(27) 



\h + k* 2 \ 

These two examples show that the phase shifts and hence the relative separation distances 
(see eq.(25)) vary as the complex parameters cc^'s change and this variation will be reflected 
in the shape of the profiles of the two interacting solitons. 



III. STATIONARY SOLITONS AND RELATIVE SEPARATION DISTANCES 

In order to realize the effect of phase shift on the shape of the solitons we consider the 
simple case of stationary limit of the two-soliton solution (5). Let us consider the situation 
in which the velocities are zero, that is kjj = 0, j = 1, 2. Further we choose a^p = a'f 1 = 0, 
a W _ e »no ? q,( 2 ) — — e V2o an( j fc nI — w here ri i0 ,i = 1,2 are real parameters (corresponding 



8 



to orthogonal modes). In this limit, the two-soliton solution(5) reduces to 
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so that each soliton is in one particular mode. 

It is also of interest to note that the stationary limit of the two-soliton solution obtained 
above is also the so-called partially coherent stationary soliton studied in the literature 
intensively in recent times [8,10] in connection with the existence of spatial solitons in 
photorefractive materials. In fact eqs. (28-32) are nothing but the stationary 2-PCS solution 
obtained in eqs. (13-15) in ref.[10]. Now we can identify the relative separation distance 
between the solitons as 
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km k 2R 2k 1R [4:kj R (k 1R + k 2R ) 



n{ki R - k 2R ) 



1 



2k 



-log 



2R 



fi{k 1R - k 2R ) 



(33) 



4k% R {ki R + k 2R )_ 

It can be very easily seen that the shape of the stationary soliton depends very much on the 
relative separation distance t\ 2 . To illustrate this in fig.(l) we have plotted (i) symmetric 
case ti2 = and (ii) asymmetric case t\ 2 ^ 0. 

One can proceed to consider the propagation of the above stationary soliton solution and 
check their shape changing properties under collision. Choosing the parameters ot^ and 
as functions of velocities (kji) such that they vanish when kji = 0, j = 1, 2, the 



nature of soliton collisions is shown in figs. 2 for the parameters = 1.0, ot[' = 1.0, 
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-2.0, fci = 1.0 + i and k 2 = 2.0 



i. 
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IV. CONCLUSION 



The soliton interactions in CNLS equations possess very rich structure. In this paper, 
we have discussed the two-soliton interaction properties of 2-CNLS equations with special 
emphasis on the nature of phase-shift encountered by solitons under collision and its depen- 
dence on the amplitudes of the modes. We have also pointed out that the much discussed 
stationary PCS solitons correspond to stationary limit of appropriate soliton solutions with 
2-PCS as an example here. Because of the complex nature of soliton interaction, multisoliton 
solution in multicomponent systems possess highly nontrivial structures. These properties 
will be presented seperately [12]. Such studies are expected to have very important appli- 
cation in optical communications, optical devices and optical computing. 
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(a) (b) 

FIG. 1: Typical stationary form of the 2-soliton solution (PCS) for the 2-CNLS system for z=0, 
see eqs. (28-32) : (a) symmetric case (t\2 = 0), (b) asymmetric case (tu = 1.0). 
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FIG. 2: Asymptotic forms of two-soliton solution (whose stationary form is similar to fig.l) of the 
integrable 2-CNLS equations (a) at z=-2 and (b) at z=2. 
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